Bubbling solutions for a planar exponential nonlinear elliptic equation with a singular source

Let $\Omega$ be a bounded domain in $\mathbb{R}^2$ with smooth boundary, we study the following elliptic Dirichlet problem $$ \begin{cases} -\Delta\upsilon= e^{\upsilon}-s\phi_1-4\pi\alpha\delta_p-h(x)\,\,\,\, \,\textrm{in}\,\,\,\,\,\Omega,\\ \upsilon=0 \quad\quad\quad\quad\quad\quad \qquad\qquad\quad\quad\,\,\,\, \textrm{on}\,\ \,\partial\Omega, \end{cases} where $s > 0$ is large parameter, $h\in C^{0,\gamma}(\overline{\Omega})$, $p\in\Omega$, $\alpha\in(-1,+\infty)\setminus\mathbb{N}$, $\delta_p$ denotes Dirac measure supported at point $p$ and $\phi_1$ positive first eigenfunction of $-\Delta\phi=\lambda\phi$ under boundary condition $\Omega$. If strict local maximum $\phi_1$, show that such has family solutions $\upsilon_s$ arbitrary $m$ bubbles accumulating to $p$, quantity \int_{\Omega}e^{\upsilon_s} \rightarrow8\pi(m+1+\alpha)\phi_1(p) \ \text{ as $s\rightarrow+\infty$.}